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Abstract 

We present a variety of series representations of the Stieltjes and related 
constants, the Stieltjes constants being the coefficients of the Laurent expansion 
of the Hurwitz zeta function £(s, a) about s = 1. Additionally we obtain series 
and integral representations of a sum S^(n) formed as an alternating binomial 
series from the Stieltjes constants. The slowly varying sum 5 7 (n) + n is an 
important subsum in application of the Li criterion for the Riemann hypothesis. 
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Introduction 

The Hurwitz zeta function, defined by £(s, a) = J2%Lo(. n + a )~ s f° r Re s > 1 and 
Re a > 0, extends to a meromorphic function in the entire complex s-plane. This 
analytic continuation to C has a simple pole of residue one. This is reflected in the 
Laurent expansion 

1 00 ( — 1 \n 

C{s, a) = + £ K —^ ln {a){s - IT, (1) 

wherein jk( a ) are designated the Stieltjes constants [2j [3l El El El [22] and 70(a) = 
—if)(a), where if) = T'/T is the digamma function. The constants 7^(a) can be written 
in the form 



( N 1 

7 fc (a) = lim V ■ — ln fc (m + a) 

N ^°° V^o m + a 



ln* +1 (JV + a)' 



N ^°° \m^o m + a k + 1 J 

In certain special cases including a > 1 a positive integer and a = 1/2 the Hurwitz 
zeta function may be expressed in terms of the Riemann zeta function ((s). In the 
case of a = 1, one puts 7^(1) = jk by convention. 

In this paper we present an array of series representations of the Stieltjes and 
related constants. Moreover we obtain series and integral representations of the sum 

4^(1)7*-, (3) 

Putting P\(t) — Bx(t — [t]) — t — [t] — 1/2 the first periodized Bernoulli polynomial 
(e.g., [111 E]) we may also write 

/oo 1 
-L^QnVdPtf), (4) 



where is the Laguerre polynomial of degree n. We have shown elsewhere jl] that 
Sry{n) is a component sum of the sum S^n) arising in application of the Li criterion for 
the Riemann hypothesis [21 E]. We have demonstrated [I] that S^(n) + n = 0{n 1 / 4 ). 
Whereas the Stieltjes constants are very closely related to derivatives of the zeta 
function at s — 1, we develop a series representation for them in terms of zeta 
function derivatives at integers j > 2. 

For a collection of formulae involving series of zeta function values we mention 
Ref. [El. 

We introduce the constants 

, / \ f°° ln j_1 (t + a) , f\t]+a + l\ , . „ 

Jo (t + a) U + a + 1 / 



Ci = r^lnfl±i > )dt=^(l). J>1- (6) 



and 

' ' ~ /l / V. t - I 
The approximate values of the first few of the latter constants are c\ ~ —0.334, 

C2 — —0.433, and C3 ~ —0.93. We have in general from the definition (6) that Cj < 0. 

We then have 

Proposition 1. For integers j > 1 we have (a) 

7i = E ^-C W H - ( - 2^)C(J + 1) - j Cj , (7) 



m=2 m 



(b) 



7i (o) = E — - J'E^ 1 )^ ln W a Li, +2 (-a) 



m=2 m fcO 



—jdj(a) — ln J+1 a 



(c) for integers n > 2 



5 7 (n) = l-ln2-(21n2 + 7-l)n + ^ 



00 Pi(x) 



X 







In f H 


--)} 


.x \ 


x) . 



roc y x 

~Ji xHl 



Pi(x) 



n 



- L^ilnx)] dx+Yl 



+ flLl_ 1 (-x)-Ll_ l{ -x/2)] J ^ )+ l 
and (d) with /(*) = \n(t + l)/t(t + 1) 

7T 2 I 



L 2 n _ 2 {^ix)dx 
(-l) m fm-2\ n 

dt, (9) 



m \m — 1 
In 



m=3 

°°Ll_ 2 (lnt) ,{[t] + r 



t 



t + 1 



~ (-l)J-i 

^ 7 1 

j= i J- 



TT 1 7 f°° 

(j'9) = - + -ln 2 2- i ln2 + ^ P,{t)f\t)dt, 



(10a) 



while with g(i, a) = ln(t + a)/(t + a — l)(t + a) 



(-l)i- 1 



00 / 1\ 

U d j) = E + "J ln ( a + 1) - ^ 



fe=i 



a 



C(2) + iln 2 (a+l) + Li 2 (-a) + i^±il 
2 2 a(a + 1) 



1 + 



1 \ f°° 

- ln(a+l)-lna+ / P 1 (t)g / (t,a)dt. 

a J Ji 

(106) 



In Eq. (8) LL,- is the polylogarithm function and in Eq. (9) 7 is the Euler constant. 

After the proof of this Proposition we present additional summatory results and 
a generalization of them. These other Propositions pertain to 70(a) and sums of 
Hurwitz zeta function values. The final section and Proposition treat the coefficients 
of the Laurent expansion of the logarithmic derivative of the Riemann zeta function 
about s — 1. 
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Proof of Proposition 1 

For part (a) we begin by forming 

co / i \m oo co / -i \m 



m=2 " b k=2m=2 



ln j k 



k=2 



( x W k x ) 
U=2 K k=2 ) 

We then apply summation by parts to write 

V [In k - ln(Jfe + 1)] ln j k = - ln([x] + 1) ln j x + j [* ln([t] + 1) 



(11) 



ln([x] + 1) ln J x + jU In(i + l)—j—dt + cf> J , (12) 



^f^ ln (M±i K ,>,. (13) 



where 



/i t \ t + 1 t 

The integral on the right side of Eq. (12) is evaluated in terms of polylogarithms 
Lij(z) = Y^k=\ zk Iki for |z| < 1 and otherwise analytically continued, by making use 
of the following. 

Lemma 1. We have 

In- 2-1 / fn — 1 



Jf ln(m)— p- eft = EC" 1 ) W e jln^~ 1 xLi, +2 (-x) + (-iy(j-l)!(l-2^)C(j + l). 

(14) 

To prove the Lemma, we make use of the representation [9] (p. 568) 

/ ln r xln(l-6x)— = (-l) r - 1 r(r + l)Li r+2 (6), Re r > -1, (15) 
J o a; 



where T is the Gamma function, and the relation [16] 

Li, +1 (-l) = -(l-2^')C(j + l). (16) 

We write 

/ ln(t + l) dt= / ln(t + l) dt- ln(t + l) dt 

Ji t Jo t Jo t 

= ln(l + a;M;)(lna; + lnM;) J - 1 (-lV(j - 1)!LL +1 (-1). (17) 

Jo w 

We then binomially expand the integrand in this equation and apply both Eqs. (15) 
and (16) and the Lemma follows. 

From Lemma 1 we may now write for Eq. (12) 



£pnJfe - ln(k + 1)] W k = - \n([x) + 1) W 



x 



k=2 



j ~ 1 (<i — 1 

-1/71 (J 1 



, ln^ 1 xLi, +2 (-x) + (-l)^!(l-2^)C(j + l)+J9(x). (18) 



From the relation between Li n (— x) and Li n (— 1/x) (e.g., [16], p. 116) 

1\ 1,„ J^ 1 ln n - 2fc x 



U n (-x) + (-l) n U n (--)=-- ln"x + 2 £ ——LiaC-l), (19) 



we may read off the asymptotic form 

Li n (-x) -ln"x + 2^ - __Li 2fc (-l), x^oo. (20) 

n\ fc=1 (n lK)\ 

Then inserting this expression into the sum on the right side of Eq. (17) gives 

/E(-l)^( j ~ X ) In^^x Li £+2 (-x) 



-/eW- 1 



ln J+1 



«+2]/2 j n <?+2-2fe 



+ 2 E 



(£ + !)(£ + 2) ^ (£ + 2-2A;)! 



Li 2fc (-1) 



J + l 



W +1 



x 



[e+2}/2 



ln- 2fc xLi 2jfc (-l). (21) 



For a given fc, the sum over £ on the right side of Eq. (21) gives zero as 



E(-i)' 

*=0 



-i ft - ^ 



£ y (£ + 2 - 2A;)! (j - 2k + 1) T(2 - 2ife) 
1 



r(2fc- l)sin7r(2ife- 1). 



(j - 2A; + 1)tt 
We now have from Eqs. (11), (12), and (21) that 

• 1 r s .. in ? /.- 



E 

m=2 



l)^)(m)= Mm E 



a;— >oo I - 1 — ' h 

k=2 ^ 



ln i+1 x + O 



ln J x^ 



x 



3 + 



- w +1 x+j Cj (x) + (-iy 3 \(i - 2-*)cu + 1)| 



lim 

t. — »r>r> \ — ' 



\n j k \n j+1 



x 



\k=2 h 



jc j + (-iyj\(i-2^)c(j + i) 



J + 1 . 

■ .in I i- lK/Iil 2 Mu./ • 1). 
In the last step we used relation (2). For the next to last step, since we have 



(22) 



(23) 



(24) 



where {£} denotes the fractional part of t, the constants Cj are certainly well defined 
and lim^oo cf^ = Cj. We have demonstrated part (a). 
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For part (b) we proceed similarly to part (a). We first have 



£ i-J2(-l)^'>(m,a) = lim if] + ~ + EM fc + a ) ~ ln ( fc + <*+!)] ^(A: + a) 



m=2 



m 



(25) 



k=0 



We next apply summation by parts, writing 

X 

^2[\n(k+a)-\n(k+a+l)} \n j (k+a) = [lna-ln(a+l)] ln J ' a+[ln(a+l)-ln([x]+a+l)] ln j (x+a) 

W-\t + a) 

— LLl{[L\ i- u -t- ±jj 

The last term here can be written 



-j / [ln(a + l)-ln([t]+a + l)]- 

Jo t + a 



-eft. 



(26) 



j f [m(a + 1) - ln([t] + a + 1)1 — ^ t + a) dt = ln(a + 1)[W a - ln j (x + a)] 
Jo t + a 



(27) 



where we have put 



.l^-( t + B ) / M+ q + l\ 
(t + a) U + a + 1 ) 



(28) 



The integral of Eq. (27) is obtained from Lemma 1 as 

ln j -\t + a) 



fX 

/ ln(t + a + l)- 
Jo 



(t + a) 



-dt 



, -/) ln ('" + 1 



W" 1 u 



ln J '- /_1 (x + a) Li £+2 (-x - a) - ln^" 1 a Li £+2 (-a)J . (29) 
We combine Eqs. (26), (27) and (29), giving 



E[hi(/c+a)-ln(/c+a+l)] ln j (k+a) = [lna-ln(a+l)] ln J a+[ln(a+l)-ln([x]+a+l)] \n j (x+a) 

k=0 



+ ln(a + l)[ln J a - \n J (x + a)} + jdf\a) 

+J EC- 1 )'" 1 ^ ( J ^ ^ t 111 '"'" 1 ^ + a ) Li ^(-^ - °) ~ ln '~"~ 1 a Li £+2 (-a)] . (30) 
We then make use of the asymptotic relation (20), obtaining 

2-_(-lKW(m,a)=hm|EnfcT^- + l11 a " ln ( " + a) + °(^T^ 



+ -A^ +1 (a; + a) + jdf (a) + J £\- 1)^! f J ^ ln^ 1 a Li, +2 (-a 

= 7j(a) + ln J+1 a + jdj(a) + j ^(— 1) £ £! ^ ~ ^ ln^" 1 a Li £+2 (-a). (31) 

Part (b) has been demonstrated. 

Remark. When a — 1, the only term remaining in the sum on the right side of 
Eq. (8) is £ — j — 1. By relation (16), we then see the proper reduction of part (b) 
to part (a). 

For the proof of part (c) we make repeated use of 
Lemma 2. We have 

gSG)^ =( - ir " 1L "- w - (32) 

that follows from the power series form of the Laguerre polynomials. We also apply 
Lemma 3. We have for integers n > 1 

c""m = i -ir n f^i n ^ xdx -(-ir s f^^ X d, + ±^. (33) 



Lemma 3 follows from induction on the integral representation [19J (p. 14) 

as) = I + A - s r ^dx. (34) 



2 s - 1 Ji x s+1 

We combine the definition (3) of S 7 (n), Lemma 2, the definition (6) of Cj, and the ex- 
pression (7) for 7j. For the ((j) values in Eq. (7) we apply the integral representation 
for Re s > 1 and Re a > 



r(s)C(s,a) 



+a— lg— (a— l)t 

o e* — 1 



-dfc. 



(35) 



at a = 1. We then obtain 



S^n) = -7n - 



n 



oo I -i \m 



^(fc-l)!V*, 



^ LlL C (fe-i)( m ) 



m=2 



m 



dx 



+ 



(e x - 1) Ji 



ooLt 2 (lnt) ln /^M + 1 



t + 1 



tZt. (36) 



We use Lemma 3 to re-express the double series on the right side of Eq. (36): 



E 

k=2 



n 



(k — 1)! \k y 



oo / i\m oo /_i\m n / \ 

E ^'H = E ^EC-i)*" 1 ? 



m=2 



m 



m=2 m fc=2 



(fe — 2)! 7i x m+1 



' ln fe 2 x dx 



(m — l) fc (A; — 1)! Ji x 



' Mn fe 1 x dx 



^m+l 



E 

m=2 



m 



°° r2 



1 X' 



rj-L n _ 2 (lnx)dx + 1 



fra- 2 
(m — 1) \m — 1 



yoo ^(3;) 



(37) 



We then make use of the elementary series 

-n E 



00 / 1 \m 

^L_LL = i_ ln2 , 



m=2 



m 



m=2 



m(m — 1) 



n — 2n In 2, 



(38a) 



10 



V (-I)" 1 _ 1 1,/, i 1\ y (~l) m _ 1 f oo. 



m=2 



m x 1 ' 



and part (c) of Proposition 1 follows. 

For part (d) we may first form using Eq. (6) 

^ roc / [i] + 1 



E ^ 0' c i) = jT ^ ln (iTf J Re « < 0. (39) 



We may note in passing that this equation properly recovers the expression for c\ at 
u — 0. At tt = — 1 we have 

N _ ^ /- fe+1 ln(/c + 1) ^ /^ln(t + l) 
An integration by parts gives the integral 



™(-iy~\ s ~ /•*+! ln(/c + 1) , /^ln(t + l), 



J*Wt_^ dt = ]nx _ ^i + I^i n ( a; + i) + 21n2. (41) 



Therefore we have from Eq. (40) 

We next apply Euler-Maclaurin summation to the sum on the right side of Eq. (36): 

lwl = r /(t)4+ 5' (i)+ r w(t)4 

= - + -ln 2 2 + i ln2 + y i P 1 (t)f(t)dt. (43) 

Combining Eqs. (42) and (43) gives Eq. (10a). 
Very similarly, we have 

^ „.j-i roo . / [t] + a + 1 \ 
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giving 

j=l k=l 

Applying Euler-Maclaurin summation then gives Eq. (10b), completing part (d) and 



E M— 0'^) = E s(*> «)-(! + - + 1) - In a. 



the Proposition. 

Remarks, (i) Concerning approximate numerical values in part (d) we have 
00 In (A; + 1) r°° 

E TT7 in - 1-25774688694, / P^t) f (t)dt ~ 0.0233638, (44a) 

fc=i ^(^ + 1) "' 1 

and 

00 (-iV'- 1 

E - — 1 0' c j) - -0-12816. (446) 

j=l J- 

(ii) The interesting sum of Eq. (44a) has a multitude of equivalent expressions, 
including a great variety of series and integral representations. By splitting the sum 
over even and odd contributions we have 



" ln(fc + l) 2 1 



\n(k + 1) ln(2A; + 1) 



fc + 1 fc 



(45) 



^1 *(* + !) ' " 2 fe tl(2fc + l) 
By using an integral representation for \nz, we have derived the representations 

~ln(/c + l) (e-* - 1) , Al f 1 (1 - x) ln(l - x) , , 

Etti tt = / ^ln(l-e"*)dt= / 46a 

^ k(k + l) Jo t y ' Jo x\nx y ' 

roo roc 

= e - *[l + ln(l - e -*)]lnt dt = -7+ / e-*lnUn(l-e-*)dt, (466) 
Jo ' Jo 

^,, + J\-. hMn{ l- e -- }dt -ll^01, (46c) 

where the second line follows from the first by integration by parts and T(x, y) is the 
incomplete Gamma function. By using the Maclaurin series for ln(l — z) in either 
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integral of Eq. (46a) we find 



h *(* + !) Sir V '• 
In addition, we have by partial summation 

"m(n + 2) ln(n + l) 



(46d) 



- Im> + 1) __ - 



n + 2 



n + 1 



(46e) 



where if n = Sfe=i 1/^ are the harmonic numbers, and we next show how this expres- 
sion may be alternatively obtained from the integral representation (46b). We use 
the generating function 



ln(l-z) = (l-z)£if n z" 



(47) 



at z — e~*. We substitute for ln(l — e~*) on the right side of Eq. (46b), interchange 
summation and integration, justified by the absolute convergence of the integral, and 
find that 

r°° °° roo 

/ e~ t \nt\n(l-e- t )dt = -y H n [ e -(™ +1 )' - e -(™+ 2 )'] \ nt dt (48a) 
Jo n=1 Jo 

\n(n + 1) ln(ra + 2 



n=l 



7 



+ 



(486) 



,n + l n + 2/ n+1 n + 2 

We note that simply if n _i — H n — 1/n and recover Eq. (46e). 

Using the generating function (47) and the right-most integral of Eq. (46a) we 
find that 

00 \ n (k A. I) 00 rl (I _ ^2 n-1 oo 

Efc^iy = E^n/ o 1 * J g dx=E^n[lnn-21n(n+l)+ln(n+2)]. (49) 



n=\ 
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The integral of Eq. (49) may be found from the special case of the Beta function 
B(x, y) integral 

/ (1 - x) 2 x q ~ l dx = B(3, q) = —. % = — + — !— , (50) 

Jo V ^ V ; ?(? + !)(? + 2) q q + l g + 2' V ; 

by integrating over q. Similarly we may determine Eq. (43) by using the other integral 
ofEq. (46a). 

Other forms of the sum of Eq. (44a) may be obtained from Eq. (46b) and the 
generating function 



oc 



9j (z) = J2n j (H n -l)z n , \z\<l. (51) 



n=l 



We may advance to gj + i from gj with the use of the operator 6 = z(d/dz): gj + x(z) 
9gj(z). We have gi(z) = — z\n(l — z)/(l — z) 2 and obtain 



£^TT)^S ,!(ff "- 1) 



In n ln(n + 1) ln(n + 2) 



(52) 



n n+1 n + 2 

One could pose the question whether all such expressions obtainable from the use of 
gj(z) are also obtainable through partial summation. 

Finally, we illustrate the application of a Binet formula for the digamma function 
[9 J (p. 943), that we use as a representation of ln(fc + 1). We have 

~ ln(fc + 1) _ ~ jjjk + 1) 1 ~ 1 ~ t dt 

hi k (k + 1) ~ hi k (k + 1) + 2 h Hk + I) 2 + hi' i f2 + ( k + !) 2 ] ( e2nt - !)" 

(53) 

We evaluate the ip sum of this equation by inserting an integral representation for it 
and interchanging operations. For the sum on the right side of Eq. (53) we apply the 
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partial fraction decomposition for the hyperbolic cotangent function [9] (p. 36). We 
obtain 



~ ln(fc + i; 



7T 



1-7 + - + 2 

12 Jo 

3 7T 2 

= - - 2 7 + — + 



1 



7T _ 1 

— coth irt - 

2t 2t 2 t 2 + 1 



tdt 



o2irt 



71 COth Tit 

t 



dt 



2 1 12 Jo 
For the last step we have used [9] (p. 328). 

(iii) The alternating sum on the right side of Eq. (9), 



3 2vrt 



-1) 



(54a) 
(546) 



00 (— 1 



m 



m-2\ n 
m — 1 



E 



-n/ (m— 1) 



(55) 



m=3 m=3 

is exponentially decreasing with n. The last integral term on the right side of Eq. (9) 
is over estimated by n — 1. 

A Corollary and further remarks 

The method of Proposition 1(b) also has relevance to Dirichlet L functions, as 
these functions may be written as a linear combination of Hurwitz zeta functions. 
For instance, for x a character modulo m and Re s > 1 we have 



k s 



m 



1 m / h 

Ex(*)( « ' 



m 



(56) 



k=l k=l 

Conversely, this relation may be inverted, yielding the Hurwitz zeta function ((s,a) 
at rational values of a in terms of a linear combination of L functions. When \ is a 
non principal character in Eq. (56) convergence obtains for Re s > 0. 
Now a suitable definition for generalized Stieltjes constants is (cf. 

ln fc+1 N \ 



lim (£^ln fc n 



m(k + 1) 



(57) 
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We first note the following from the first equality in (56). We have 

oo / i \n 00 00 / i\n 

E [ —^hiyL«\n, x) = E X(k) ln J k £ ^ 



71=2 



E *(*) ln j k 



k=2 



k=2 



1 . (k+l 

— in 



n=2 



fc 



fc 



Putting 



= Jim ^ x(k)—— + E *(*0 t ln k ~ ln (^ + !)] lnJ fc 

U=2 ^ k=2 



A(x) = J2x(k)[^k-Hk + l)], 

k=2 



(58) 



(59) 



and applying summation by parts we have 



V X (k) [In k - \n(k + 1)] ln J k = A(x) \n j x-j [* A(t) 

k=2 ' Jl 



In'" 1 1 
t 



dt. 



(60) 



Alternatively, we now proceed in terms of Hurwitz zeta functions and write from 
the product rule and Eq. (56) 



(-1)^)( S)X )=E(-1/ 



'j\ ln J m 



E x(k)c (e) k £- ] ■ 



A: 



fc=i 



m 



(61) 



We then have 
-1) 



1) 



n J 



— 1 n n=2 n «-« 



n=2 



£=0 



^EX(A:)C W 



n, 



fe=i 



m , 



-\) n * (j\\n j - e m 



n=2 



£=0 



ExWE 



00 J /A W^m m 



EE 

9=0 ^=0 



k=l q=0 
1 



k=i 



mq + k 



\n £ (q + fc/m) 
(g + k/m) n 

mq + k + 1 
mg + A: 



In 



(62) 
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The sum on q here is 



In (q + k/m) 

q=0 



1 



mq + k 



-In 



mq + k + 1 
mq + k 



. £ In (q + fc/m) + + q _ ln ( mg + jfe + !)] \ n \ q + fc/ m ) [, . ( 63 ) 



q ^ mq + k q=0 
Summing by parts, we have 



^[ln(m<7 + k) - ln(mg + k + l)) ln £ (q + k/m) = [In A; — ln(k + 1)] \n e (k/ 



m) 



q=0 



-B(x) \n e (x + k/m) - if* B(t) ] - 

Jo 



ln' -1 (t + A;/m) 
t + k/m 



dt. 



(64) 



where we have put 



B(x) = y^[ln(mg + k) — ln(mg + k + 1)]. 

9=1 



(65) 



The integral term of Eq. (64) gives an extension of the integral of Eq. (26) when 
m > 1. 

From Proposition 1(b) we obtain 
Corollary 1. We have for Re a > and Re b > 



In 



r(6) 



oo 

+ ln(A;+a 



POO 

b In 2 b-a In 2 a+ / [l+ln(a+t)l In 
Jo 

k + a+1 



'[t)+a + l 



k=0 



1 - (k + a) In 



k + a 



^ t + a + 

oo 

H k + b ) 



i \ f°° 

Yj dt -J Q [l+ln(fe+0]ln 



"M+6+T 
£ + 6 + 1 



fc=0 



1 - (Jfe + 6) In 



k + b 



+6+ln6[-6+(H-6)ln(6+l)]+Li 2 (-6)-a+lno[o-(H-a)ln(o+l)]-Li 2 (-a). (66) 
In contrast, a standard sum representation is [H] (p. 936) 



In 



r(o) 



In 



+ E 

fc=l 



(6 - a) In 1 



In 



+ a 



(67) 
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Proof. We use the summatory relation [8] (Proposition 1) 

T(6) 



oo y 

-ihn+i(a) - 7«+i(&)] = In 



n=0 



r(o) 



Re a > 0, Re b > 0. 



(68) 



We write 



m=2 



:-iy 



C (i) (m, a) - C ij) (m, &)] = £ ln J (fc + a) 

fc=0 



1 , /fc + a + 1" 
— m 



k + a 



k + a 



fc=0 



In 



k + b + l 



(69) 



k + 6 \ fc + 6 

Putting j — n + 1 and summing Y^^=o(^/ n ^) gives the terms of line 2 of the Corollary. 
The integral terms of the Corollary result from using the definition (5) of d n +i(y). 
For the polylogarithm terms we have 

£^±(-W")w^«um(-«o 



n=0 



€=0 



£(-l)^!Li, +2 (-a) E 



(n + 1) " /n 



£=0 



77 ' 



n-£ 



aE(-l)'(l + lna + £)Li, +2 (-a) 

oo 



1 In a 

+ 



(j + 1) 2 J(J + 1) 



-a + lna[— a + (1 + a) ln(a + 1)] + Li2(— a). 



(70) 



Combining terms completes the Corollary. 

Akin to Eq. (24) we have from Eqs. (5) and (6) the weak upper asymptotic 
bounds 

\ Cj \ = 0[(j-l)\], \d j (a)\ = 0\T(j,]na)], j - oo. (71) 



These bounds miss cancellation, so are large over estimates. 

Finally in closing this section we write a summation form for the constants dj(a). 
From Eq. (5) and the use of Lemma 1 we have 

rk+a+l in^ 1 11 

dj(a) = 22 [\n(k + a + 1) - \n(u + 1)] du 

i, nJ k+a U 



k=0 



1 oo 

- [W(k + a + 1) - ln J (A; + a)] ln(k + a + 1) 

3 k=o 



fc=0 €=0 



hi 3 - e -\k + a + l)Li^ +2 (-fc - a - 1) - ln^" 1 ^ + a)U e+2 (-k - a) 

(72) 



Other expressions for 70(a) 



Recalling that 70(a) = —if) (a), we have 
Proposition 2. We have (a) 



70(0) = — + ln(o + 1) - £ { 



m=2 



m 



C(m,a) - 



(b) 



00 ( 1 \m 



m=2 

and for \z\ < |a| (c) 



rC( m ) a) = —a In a + In T(a + 1) — -ip(a) — - ln(27r) + a, 

771+1 2 2 



(73) 



(74) 



00 ^rn r 



E 

m=2 



rn 



C(m,a) 



1 



- + In (l - - J - lnT(a - z) + lnT(a) - z^(a). (75) 



Proof. For part (a), we form 

:-i) 



00 l 1 \m 



m=2 



C(m, a) 



00 / -1 \m oo 



m=2 
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= E 

k=l 



k + a 
ln(l + a) + lim 



ln(/c + a + 1) + ln(fc + a) 



E 



-ln([x] + a + l) 



_ fc=1 k + a 

= ln(l + a) - ^(a + 1) + Jim [ip([x] + a + 1) - ln([ar] + a + 1)] 
= ln(l + a) -^(o+l)+|mO ' 
= ln(l + a) - ^(a) 



(76) 



z->oc y [x] + a + 1 ^ 
1 

a 



(77) 



In the last step we have used the functional equation of the digamma function. 

An alternative method of proving part (a) is by applying the integral representa- 
tion (35). Then we have 



E 

m=2 



{-iy 



m 



C(m,a) - 



1 



E 

m=2 



-iy 



m 



ll/o 



oo ^m-lg-(o-l)t ]_ 

dt 



(in — 1) ! jo e* — 1 



1 , / IN ~ (-1) 
- + ln 1 + - -E 



m=2 



to! Jo 



^m— lg— (a— l)t 



= -i + ln(l + i 

a V a 



oo g-(o-l)t Y 

- r - j --(e- t + *-l)dt 
o e 1 — 1 r 



— + In (l + - In a + ^>(a) 
a V a/ 

= -- + ln(a + 1) + ip(a). 
a 



(78) 



The absolute convergence of the integral representation (35) justifies the interchange 
of summation and integration above. 
For part (b), we have 

oo / i\ra oo / T\m oo i 

lib T 1 ___ 9 



m=2 



f 2 (m + 1) f (* + a)' 
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00 [ 111 
= 2|(fc + a)pn(Ar + o+l)-ln(* + a)] + 2(jfc^- 1 j 

lim ^{(k + a + l) \n(k + a + 1) - (k + a) ln(A; + a) - ln(A; + a + 1) + 



1 1 



2 (jfe + a) 
(79) 



lim {([ar] + a + l)ln([x] + a + 1) - alna + lnr(a + 1) -lnTQa] + a + 2) 



+ _^([ x ] +a + l)_^( G )]_([ x] + l)j. 

We then apply Stirling's formula for lnr(z) and part (b) follows. 

For part (c), we proceed as in the second proof of part (a). We have 



(80) 



oo m 

T.- 



m=2 



m 



C(m, a) 

v ' a m . 



oo m 

— 



m=2 



m 



(m — 1)! Jo e* — 1 



^m-lg— (a— l)t j 

di 



z / z \ r°° e _(a_1) * 1 

- + In 1 / — r— I 

a \ aj Jo (e* — 1) t 



e 2 * - zt - l)dt. 



We use p] (p. 939) 

lnT(a — z) — lnT(a) 



^ — {a — z)t g— at 

1 - e-* 



-t 



tit 

7' 



52) 



along with an integral representation of 7 plus the digamma function [9] (p. 943) and 
of the Euler constant [9J (p. 946) in the form 



7 



,-t 



1-e- 



dt. 



(83) 



Proposition 2 is finished. 

Remark. Part (a) of this Proposition gives the 1 extension of Ref . [20] . Part 
(b) gives the a 7^ 1 extension of Refs. [15] and [TB] . 
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Generalizations of sum of Hurwitz zeta function values 

Here we generalize the sums of the previous section, showing 

Proposition 3. For integers j > 1 we have 

00 (—i) m x r / 1 

E ^TTCK a) = (-iy Mm £ { (* + a)''" 1 " (* + «)'' ln 1 + r~T 

m =2 J fc=0 L \ At -t- a 

- (H *I,^(H<)' <84) 
The proof of this Proposition depends upon 

Lemma 4. For integers j > and k e C let 

M,-(ib) = V ; ' r- — . (85) 
Then we have the expressions 

= (7T2) (F^F 2Fl + 2 ^ + 3! (86) 

= (-!)>(* + ayM (k) - (-l) J (k + aY £ 1 J , (87) 

m =2 m l, K "+" a ) 

where 2 F\ is the Gauss hypergeometric function and 

M »w = ^- ln ( 1 + ^)' (88) 

For the proof of Eq. (86), let (a) n = T(a + n)/T(a) be the Pochhammer symbol. 
By noting that (j + 2) m / (j + 3) m = (j + 2)/(j + m + 2) we have 



m 



- (m + j + 2) (A; + a) m + 2 (j + 2)(k + a) 2 ^ m! (j + 3) m V fc + a 

(89) 
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By the series definition of the function 2 F±, Eq. (86) holds. In particular, M (k) is a 
known special case. 

For Eq. (87), we observe by shifting the summation index that 

M j+1 (k) = —(k + a)M,(A;) + ^L- (90) 

In addition to this recursion relation we have 

^(^(-m^ay^t^——, ( 91 ) 

and Eq. (87) follows. 

For the proof of Proposition 3, we write 

oo / -r \m oo 

E L ^-CKa) = E^W ? (92) 

m=2 m ~r J k=0 

and then apply Lemma 4. 

The manner in which the sums of terms rational in k + a in Proposition 3 may be 
evaluated in terms of polygamma functions ip^ is given by 
Lemma 5. We have (a) 

x I i \m—j 

J2(k + ay- m = ^-M — U (m " J_1) («) - ^ (m_j_1) (N + a + 1) , (93) 

fe =o (m-j - 1)! L 1 

and (b) the asymptotic value of the second term on the right side of Eq. (93) is given 
by 

(-l) m ~ J ,(^-1)1, , ^ 1 1 Q ( 1 \ 

(m-j-l)r Wnui-i) ( m - j - i) ([ x ] + a + \([x] + a + l) m ~i ) ' 

x — > cx). (94) 
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For the proof of part (a) of this Lemma, we note that 

J2(k + ay~ m = C(m - j, a) - C(m - j, [x]+a + l). (95) 

k=0 

We then use the relation 

(_l)n+l 

C(n + l,x) = ^j ^ {n) (x), n>l. (96) 

For part (b) we apply the known asymptotic form of the polygamma functions [1] . 

Remarks. The asymptotic evaluation of terms of the form (k + a) 3 ' [In (A; + a + 1) — 
ln(/c + a)] in Proposition 3 may be performed by way of Euler-Maclaurin summation. 
So at j = 1 one way to write the result is in terms of the Glaisher-Kinkelin constant. 
Similarly, at j = 2 and j = 3 the results may be written in terms of constants B and 
C that have previously appeared in the literature (e.g., [16], pp. 36 and 37). 

Examination of Eq. (95) shows that j there need not be an integer. We may 
therefore differentiate it with respect to j and find 

X 

J2(k + a) j - m \n(k + a)= ( f (m - j, [x] + a + 1) - (\m - j, a), (97) 

k=0 

where ' denotes differentiation with respect to the first argument. Also for use in 
connection with Proposition 3 we have by the integral representation (35) and Euler- 
Maclaurin summation ((s, a + x) = x x ~ s j (s — 1) + 0(x~ s ) as x — * oo, giving 

('(s,a + x) = lnx — — + 0(x~ s lnx), x — > oo. (98) 

s — 1 (s — \y 

Additionally we may further generalize Proposition 3 to 
Proposition 4. We have for Re a > 0, j > —1, and \t\ < 1 (a) 

]T(-l) fc C(fc,a)7-p- = I x j [^(x + a)-ij(a)]dx, (99) 
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and (b) 

t k+j 



~f~ 'J ft 

(-l) m J2(-l) k (k) m ((k + m,a)-—= / x j [^ m \x + a)-^ m \a)]dx, (100) 
k=2 k + j Jo 

or, equivalent ly, (c) 

V(-l) fc (fc) m C(A; + m,a)- =m! / x J [C(m + 1, x + a) - ((m + 1, a)]cfe. (101) 

fc=2 k + j Jo 

Part (c) follows from part (b) by the use of relation (96). For the proof of part 
(a) we start with the identity 

OO 

£(-l) fc C(*:, a)^ 1 = 1>(x + a) - ^(a). (102) 

k=2 

We then multiply both sides of this equation by re- 7 and integrate. The condition on 
j is imposed due to the fact that tp(x) = — 1/x — 7 + 0(x) as x — > 0. For part (b) we 
repeatedly differentiate Eq. (75) with respect to a, use 

^C(k, a) = (-l) m (k) m ak + m,a), (103) 

and the Proposition is completed. 

Remark. Proposition 4 may be extended by multiplying Eq. (102) by other 
functions and integrating over other intervals. For instance we have 
00 r-t r-t 

y2(-l) k ((k,a) / x k+ P-\l-xy i {l-zx)- a dx = / xP{l-xy{l-zx)- a [4){x+a)-i){a)}dx. 

(104) 

Once again such a result may be differentiated with respect to a and Eqs. (103) and 
(96) used in order to obtain certain integrals of differences of Hurwitz zeta functions. 
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Moreover, Eq. (100) could again be multiplied by say t n and integrated. For special 
values of t including 1 explicit results are obtainable by using various representations 
of the digamma function. 

By using various integral representations of the Hurwitz zeta function we may 
determine evaluations for the integrals of Proposition 4. We illustrate this in relation 
to Eq. (92). If we put for integers j > 0, k E C, and \t\ < 1, 

oo ( '\\ m £ m 



MAk, t,a) = . 



(105) 



then we have 



Mj(k,t,a) 



( 3 + 2) (k + af 

iyr j {k + a y 



2 F 1 (l,j + 2;j + 3; 1 



k + a 



i+l / i \ mj.m 



and 



M0(M ' O) = ^- ln ( 1 + ^)- 



(106) 
(107) 

(108) 



These expressions follow in like manner to Lemma 4. We then have 
Proposition 5. For Re a > 0, integers j > 0, and \t\ < 1 we have 



k=2 



f 

JO 



t 2 



1 



,2 (J + 1), 



(x + a) 2 (x + t + a) 
t 

J 



jMj(x,t,a) 



P\(x)dx 



(109) 



M j (p,t,a) + —- -ln(l + */a)-2 / ]T 



00 ^ sin(/c tan -1 y/a) dy 



( Q 2 + y 2)fc/2 ( e 2^ _ X) 

(HO) 
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= {l~ (JTT) ) Mj (0 ' *' a) + JTi ln(1 + t/a) 

r°° / 1 1 1\ e~ aw 

+ / 1 9 M J (0,t,l/w)dw. (Ill) 

Jo VI — e *° u> 2/ w 

In the proof of Eq. (109) we use the representation 

a~ s a 1 ' 3 f°° PAx) 

C(s,a) = — + --a r -±J—dx, Res>-1, 112 

2 s — 1 Jo (x + a) s+l 

that is based upon Euler-Maclaurin summation. For the proof of Eq. (93) we use 
Hermite's formula for the Hurwitz zeta function (e.g., [21J p. 270), 

a~ s a 1_s r°° sm(s tan" 1 via) dy 

Similarly, we may use a different integral representation for ((k,a) to obtain Eq. 
(111). We however describe a second method based upon the relation (96), whereby 
we obtain 

OO -4-k OO 4-1' 1 

g ( _ 1)tc(M) _ = g w . (114) 

We make use of a classical Binet formula for the digamma function [16J. Upon 
repeated differentiation, we find 



XB _!(n-l)! (-l) n n! 



^W(^) = (-i)"- 1 ^ — ^ ;- -(-iy 







1 1 



1 - e~* t 2 



t n e~ tz dt, Re z > 0. 

(115) 

Based upon this expression, inserting 

^ (fc_1) (a) into Eq. (114) and again interchanging 
summation and integration we determine Eq. (111). The interchange is justified by 
the absolute convergence of the integral of Eq. (115) and Proposition 5 is completed. 
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Remarks. It is possible to re-express the summation on the right side of Eq. (110). 
It is also possible to insert Binet's formula for ip and into Eqs. (99) and (100) 
and thereby obtain other expressions. This too we omit. 

An extension and brief discussion 

In principle, a result of the shape of Proposition 1 holds for all convergent Dirichlet 
series. Rather than describe this generally, we illustrate this for a Dirichlet series of 
particular interest. Namely, for 

(' ~ A(jfc) , . 

7( s ) = -E4lA Re S >l, (116) 

^ k=2 K 

where A is the von Mangoldt function [TTIHU] . Like ( itself, the logarithmic derivative 
(116) has a simple pole at s = 1 and the Laurent expansion analogous to Eq. (1) is 

77^ = -— T-EVpis-lf, \s-l\<3, (117) 

^V b ) b 1 p=0 

wherein the constants rjj can be written as 

[-l) k 1 \n k+1 N 



Vk 



k\ n 



r 1 * / k \n k+1 N\ 
hm > — A(m) \n k m . (118) 



The r)j coefficients are related to the Stieltjes constants 7^ by a simple recurrence 
relation (e.g., [7J, Appendix). 

We now have 
Proposition 6. For integers j > we have 



1 - -1 ™ c'\ { \ \ -1 y 

m = -- V - — -M m + ^— -L- hm 



]T A(fc) Wkln 

k=2 





ln j+1 N~ 


m 


J + 1 . 



(119) 
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and 

Corollary 2. We have 

'k + r 



770 = -7 = - £ y( m ) + i im 



£A(fc)ln — ^ — \nN 



.fc=2 



(120) 



On the right side of Eq. (119), the dominant asymptotic order ln 3 ' + N/(J + 1) 
is cancelled by the sum, leading to a new set of constants analogous to the second 
and third terms on the right side of Eq. (7), whose explicit determination remains 
open. Again partial summation may be applied to the right side of Eq. (119) and 
alternative forms written. As we do not wish to invoke the Riemann hypothesis, we 
do not pursue this here. 

An approximate value of the first contribution on the right side of Eq. (120) is 
0.2419488514703058, enabling a numerical approximation for the constant coming 
from the limit term. In comparison, the term (' (2) / [2( (2)} = (7 — In 2 + hi7r + 
12 In A) /2 ~ -0.2849804965472664 where A is Glaisher's constant. 

For the proof of Proposition 6 we form 



00 



1 g (^r / g\ {J) = j_ry_ 



:i2D 



We then add the expression (118) and Eq. (119) follows. 

The importance of a new representation of the rjj coefficients lies especially with 
new representation and estimation of the critical alternating binomial sum 

s,(n) = -f; ( n Wi = E(-i) m ( n Wil- ( 122 ) 

The linearity or sublinearity of ^(n) = Sj(n) + Sx{n) in n, wherein the oscillating 
sum S^in) involves logarithms and values of the von Mangoldt function [4], more 
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than suffices to verify the Riemann hypothesis under the Li criterion. Indeed, the 
Riemann hypothesis fails only if S\(n) and hence S^n) becomes exponentially large 
in n and negative. 
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